Rabi oscillations of a qubit coupled to a two-level system 
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The problem of Rabi oscillations in a qubit coupled to a fluctuator and in contact with a heath 
bath is considered. A scheme is developed for taking into account both phase and energy relaxation 
in a phenomenological way, while taking full account of the quantum dynamics of the four-level 
system subject to a driving AC field. Significant suppression of the Rabi oscillations is found when 
the qubit and fluctuator are close to resonance. The effect of the fluctuator state on the read-out 
signal is discussed. This effect is shown to modify the observed signal significantly. This may be 
relevant to recent experiments by Simmonds et al. [Phys. Rev. Lett. 93, 077003 (2004)]. 
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Recent experiments have demonstrated Rabi oscilla- 
tions in macroscopic quantum systems (qubits). 0,0,0, 
0000 These oscillations decay rather quickly due to 
interaction with the environment even in pure systems. 
An interesting behavior was observed in spectroscopic ex- 
periments with Josephson qubits. 0] Both the decay rate 
and the oscillation pattern were strongly dependent on 
the qubit cigenfrcquency ramped by external parameters. 
Namely the oscillations were significantly suppressed in 
the vicinity of certain eigenfrcquencies. This was inter- 
preted 0] as an influence of some two-level systems (fluc- 
tuators) located in the qubit environment. The suppres- 
sion is strong when the fluctuator's and qubit's energy 
splittings are very close. 

Below we present a theory of a qubit interacting with 
an external AC field and a single fluctuator coupled to 
the qubit. In addition, we assume that the system inter- 
acts with some thermal bath providing both phase and 
energy relaxation. The account of the energy relaxation 
distinguishes our model from a similar approach used in 
Ref. whereas in the recent numerical work of Ref. 
no account is taken of the phase relaxation. By solving 
the kinetic equation for the density matrix we compute 
level populations versus time. The results explain the 
strong influence of the resonant fluctuator on the Rabi 
oscillations of the qubit in agreement with experimental 
findings. 

We will characterize the qubit by a spin S interacting 
with a fluctuator represented by a spin s. The Hamilto- 
nian in the form 



H{t) = Hq + Hf + H q -f + H m an(i) 



(1) 



takes into account the qubit and the fluctuator (H q 
and 7if) as well as the interaction between them (Ti. q _{). 
We also included into 7i(t) the "manipulation" part 
7^man(£) = S X F cos tot, which in spin terms is an oscil- 
lating magnetic field in the x-direction applied to the 



qubit (F is the Rabi frequency). The other parts of the 
Hamiltonian can be written as 



7~tq = % 



"S Z ,H q _f = -{S X S X + SySy) .(2) 



Here E is the distance between the qubit levels, e is the 
distance between the fluctuator's levels, Si and Si are 
the Pauli matrices acting respectively in the spaces of 
the qubit and fluctuator, u is the off-diagonal coupling 
constant. ^3] Below we will use the so-called resonant 
approximation replacing cosujt — > (1/2) exp(iwt), see e.g 
SlichterHiJ. This approximation, which neglects higher 
harmonics of the response, is valid close to the resonance, 
i. e., when \E — u>\, \e — u)\, \u\ <C u>. We also omit the 
diagonal qubit-fluctuator interaction cx S z s z which is im- 
portant only when the qubit and the fluctuator are far 
from the resonance. On the contrary, the off-diagonal 
part is important only when the qubit and the fluctuator 
have close energy splitting, \E — e| <§C E, e. The impor- 
tance of this interaction was stressed in Ref. We will 
also assume that 



<^ < T, 



(3) 



where u and T are measured in energy units. Only in 
this case the qubit acts as a resonant system. 

We are not going to take the thermal bath explicitly 
into account. Instead we introduce damping phcnomeno- 
logically into the equation for the density matrix evolu- 
tion. 

The Hamiltonian in the resonant approximation 
can be expressed as a 4 x 4 matrix 



n 



Fe luJt \ 

-E + e u Fe lult 

Fe~ lut u* E-e 
Fe -iwt 



(4) 



\ Fe- Wt E + e j 
Shown in Fig.n arc the energy terms of the Hamiltonian 
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FIG. 1: Diagram of the terms of the 4-level system consist- 
ing of the qubit and the fluctuator. Shown by arrows are 
transitions induced by the AC field. 



versus the qubit energy splitting E for a fixed fluctuator 
splitting e and u = O.le. Near the resonance, E = e, all 4 
levels are involved in the AC-induced transitions. In this 
region one can expect strong influence of the fluctuator 
on the qubit response. 

The 4-level system qubit+fluctuator is characterized 
by a 4x4 density matrix Pfj, v {t) which diagonal matrix 
elements Tin, n lh n tl> n TT describe the occupations of 
each of the levels. In the presence of the high-frequency 
driving force cos uit the off-diagonal matrix elements also 
quickly oscillate: 



if*e~ iut 
ig* e~ wt 
\i)*e- 2wt ig* 



ik* 



-ig e 
—ik 



-ile iu 
-im e 1 ' 
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(5) 



It is convenient to introduce a matrix with elements that 
all vary slowly with time: 



if* fin 
ig* ik* n n 
V ij* il* im* 



f -ig -ij \ 
ik —il 



—im 



(6) 



The matrix © is just the density matrix in the frame 
rotating with frequency to around the z-axis in the qubit 
spin space. In this frame the von Neumann equation, 
dtp = i[H, p], involves a time-independent Hamiltonian, 
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Here £ = (E - w)/F, e = (e - w)/F, tj = u/F. The 
matrix elements Hv.\ = H31 and H-u = 'H'h correspond 
to transitions between the qubit levels, while H23 = W| 2 
describe simultaneous "flips" of the qubit and the fluctu- 
ator. 



Our next step is to rewrite the von Neumann equation 

p = Lp (8) 

as an equation for the vector p with 16 components 

{nn 1 n iT ■> n Ti' n TT' 9, 9*, U* f, f* m,m*,k, k*, j, j*}. (9) 

The 16x16 evolution matrix L is read off from the com- 
mutator [H, p\. 

When no heat bath is coupled to the 4-level system, the 
matrix L has rank 12. This reflects the fact that in the 
absence of dissipation there are four conserved quantities, 
namely the occupation numbers of the eigenstates of the 
Hamiltonian Eq. J7J). 

The system reaches a stationary state only in the pres- 
ence of dissipation. We take the dissipation into account 
in the following way. First, we assume that at large 
time the system will reach a stationary, however non- 
equilibrium, state p°° where the matrix elements depend 
on the relationship between the pumping amplitude, F, 
and the phase and energy relaxation rates. Given the 
"dissipation operator" T, which takes into account both 
phase and energy relaxation one can determine p°° as a 
solution of the equation 



{L + t)p°° = 0. 



(10) 



The explicit form of the dissipation operator T is de- 
termined by the particular model chosen for the thermal 
bath. At the same time, the diagonal matrix elements of 
the stationary density matrix 



Pdi; 



(11) 



can be specified phenomenologically in a model- 
independent way. The dynamics of the remaining 12 ma- 
trix elements p^g-.^iag i s rather insensitive to the details 
of the model. 

The situation here is somewhat similar to nonlinear AC 
conductivity in disordered metals: when the scattering of 
itinerant electrons is predominantly elastic, the details of 
the inelastic relaxation are not important in spite of the 
Joule heat. This inelastic relaxation rate enters usually 
only though the effective field-dependent temperature, 
while the temperature-dependent conductivity can be de- 
termined without taking energy relaxation into account. 

Adopting this strategy we choose an approximate dis- 
sipation operator (16x16 matrix) r and present the op- 
erator L + Tq as 
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The off-diagonal and diagonal elements of the density 
matrix p°° are connected through the 12x12 matrix Mo 
and the 12x4 matrix M\: 

PrfF-diag = ~ M -^lPdiag ■ (I 3 ) 

The advantage of this method is that as soon as sta- 
tionary level populations are specified the results only 
weakly depend on the specific form of the matrix To . For 
simplicity we adopt a model similar to the one conven- 
tionally used to describe relaxation in various two-level 
systems, see, e. g., Ref. If the qubit does not directly 
interact with the environment, while the fluctuator has 
a finite dephasing rate 7 and energy relaxation rate 27, 
then r is diagonal (Tq)^ = Ji5ik where 

_ f -27, 1 < i < 4 
ll ~ \ -7 , 9 < i < 12 . 

This damping matrix differs from the one used in Ref. 
by taking into account relaxation of the diagonal ele- 
ments of the density matrix. 

Given p°°, the solution of the equation p(t) = (L + 
To)p(t) can be presented as an expansion 

16 

p(t) = p x +Y,c- l e t e^ t (14) 

i=2 

over the eigenvectors, ej, of the matrix L + Tq, Xi being 
the corresponding eigenvalues. For weak dissipation 
and Xi are close correspondingly to e[ ^ and A 2 - ^ - eigen- 
vectors and eigenvalues of L. We start summation in Eq. 
(|14(l from i = 2 to emphasize that Ai = regardless of 
the dissipation if we choose e.\ such that its scalar prod- 
uct with p is trp = nn + + n^i + n^. Accordingly, 
the part of p(t) proportional to e\ should be included 
into p°°. Note that in general p°° can not be presented 
as a linear combination only of ej with i < 4. 

One can specify the initial conditions assuming that 
the external AC field is switched on at t = 0: 

Pu = n u (0) = l, ft^ wl (0)=0. (15) 

We have checked numerically that the results are only 
weakly sensitive to initial conditions. It is more impor- 
tant to determine p°°, which strongly depends on the 
pumping strength, F, as well as on the relaxation times 
of the phase, T2, and energy, t\. For strong pumping 
F 2 T\T2 3> 1 the levels with different orientations of the 
qubit have almost the same occupations. The stationary 
populations of the fluctuator level depend on the inter- 
play between its coupling with the qubit and the energy 
relaxation rate. If the coupling is weak, the stationary 
occupancies of the levels are given by a thermal distribu- 
tion: 

n TT = n iT ■> n U = n ll 1 n 1l/ n U = e~ e/kT ■ 



Since we do not focus on the temperature dependence we 
will set the ratio n-j-f /tin equal to some constant which in 
principle depends on the Rabi frequency F, interaction 
strength u and the fluctuator relaxation rates. 

We have calculated the probability to find the qubit in 
its upper state, n^(t) + n^i(t), from Eq. (|H|) and plotted 
it as a function of time in Fig. [5] 
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FIG. 2: The qubit upper level occupation, n T j. + n TT , as a 
function of time measured in Rabi periods, _F _1 . The param- 
eters are: 7 = 0.03F, 77 = 0.1, £ = and n-ff/n-fi = 0.43. 
Values of the fluctuator's deviation from the resonance, e, are 
shown near the curves. 

The qubit is assumed to be in resonance with the 
pumping field, £ = 0. We observe clear Rabi oscillations 
with frequency F and a decay in time due to finite relax- 
ation, 7 = 0.03-F 1 The oscillations are perturbed by the 
fluctuator coupled to the qubit with the strength 77 = 0.1. 
Two curves correspond to different shifts of the fluctua- 
tor from the resonance, e = ±0.3 (these two curves are 
practically the same) and e = 0.01. One can see that 
close to the resonance there are beats in the Rabi os- 
cillations. This should be expected from Fig. 2] since 
at the resonance both transitions participate. Also, for 
times shorter than the inverse beating frequency or larger 
damping the resonance effect amounts to a suppression 
of the Rabi oscillations. It is clear that the closer the 
fluctuator is to resonance the stronger it suppresses Rabi 
oscillations. This fact confirms that suppression of Rabi 
oscillations observed in Ref. Q] can indeed be induced by 
coupling to a fluctuator close to resonance. Large damp- 
ing, certainly, makes the beats less visible. 

Experimentally, the qubit state is probed spectroscopi- 
cally by inducing transitions from the upper state to some 
high-energy state having small lifetime due to tunnel- 
ing. One can start from the assumption that the tran- 
sition probability is independent of the fluctuator state. 
This is why we present the results for the sum + n^ i 
in Fig. |2 see also Ref. H. 

In reality, the transition probabilities from the states 
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FIG. 3: Energy levels and transitions. 



HI) and HI) can significantly differ. Moreover, transi- 
tions from the qubit down states |||) and |J.|) can be 
induces. For example, if the qubit is resonant with the 
fluctuator, the levels ||J.) and |||) will be close in energy, 
while the level |||) will be removed from this close pair 
(see Fig. 01. 

It is then hard to expect that the driving force is reso- 
nant with transitions from the levels |||) and |||) without 
also being coupled to the |||) state. Different response of 
the different levels to the measurement pulse can signif- 
icantly change the shape of the curve. To illustrate this 



In conclusion, we develop a general scheme of analyt- 
ical analysis of time-evolution of the density matrix of a 
multi-level quantum system coupled to a heat bath. In 
this scheme the dissipation is taken into account phe- 
nomenologically, while the quantum dynamics in the 
presence of arbitrarily strong driving force is evaluated 
explicitly. We applied this scheme to the system of a 
qubit coupled to a fluctuator which probably was exper- 
imentally realized in RefQ We demonstrated that if the 
fluctuator is close to resonance with the qubit, the Rabi 
oscillations of the qubit are suppressed at short times and 
demonstrate beatings when the damping is weak enough. 
We also found that if the read-out signal depends on the 
state of the fluctuator, the visibility of the Rabi oscilla- 
tions can be substantially reduced. This effect can natu- 
rally explain the experimental results of Ref Q. 
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FIG. 4: Solid curve - weighted average of the qubit upper state 
population, 0.5rif| + riff, for e = 0.01. Other parameters are 
the same as in Fig. 0. For comparison the sum rifj. + ^TT i s 
shown by dashed curve. 

point we plot in Fig. ^a weighted average, 0.5rif j + riff, 
assuming that the escape from the state |||) is twice as 
probable as the escape from the state || J.) (solid curve). 
Interestingly, this curve is much more similar to the ex- 
perimentally observed one Q than the simple sum shown 
for comparison (dashed line). Namely, the plotted quan- 
tity oscillates never exceeding 1/2 and tends to « 1/2 at 
large time. As far as we understand it, it is unclear what 
combination of the occupation numbers of the states is 
actually probed in the experiment. In the second set of 
experiments a different readout scheme was applied, 
but no Rabi oscillation traces where presented. 
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